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Abstract 

The conductance of a contact, having a radius smaller than the Fermi wave length, on the surface 
of a thin metal film is investigated theoretically. It is shown that quantization of the electron energy 
spectrum in the film leads to a step-like dependence of differential conductance G (V) as a function 
of applied bias eV . The distance between neighboring steps in eV equals the energy level spacing 
due to size quantization. We demonstrate that a study of G (V) for both signs of the voltage maps 
the spectrum of energy levels above and below Fermi surface in scanning tunneling experiments. 

PACS numbers: 74.55. +v, 85.30.Hi, 73.50.-h 
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I. INTRODUCTION 



Today a fairly large number of papers have addressed the problem of calculating point- 
contact conductance for use in analyzing and interpreting scanning tunneling microscopy 
(STM) experiments (for reviews see, for example, [H, 0]). The low symmetry of the problem 
and the wide variety of objects under study do not allow developing a general theory of STM, 
and different approaches for specific problems are used. The theory papers on this subject 
can be divided into two groups: One uses methods taking into account the specific atomic 
structure of the STM tip and that of the test specimen. These methods make it possible 
to reproduce the crystallographic structure of the sample surface in the calculated STM 
images and this is very useful for arriving at a correct interpretation of experimental data. 
The main deficiency of this approach is the lack of analytical formulas for the STM current- 
voltage characteristics as numerical calculations must be performed for every specific case. 
The other group of works exploit simplified models of noninteracting electrons which allows 
finding relatively simple analytical expressions that describe the STM current qualitatively. 
For this reason such theoretical results are widely used by experimentalists. 

One of the first free-electron models describing STM experiments was proposed by Tersoff 
and Hamann [3| whose theoretical analysis of tunnel current is based on Bardeen's formalism 
[ij , in which a tunneling matrix element is expressed by means of independent wave functions 
for the tip and the sample within the barrier region. Using the model wave functions the 
authors [3] showed that the conductance of the system is proportional to the local density 
of states of the sample at the tip position. In principle it is possible to extract information 
on subsurface objects (single defects, clusters, interfaces etc.) by STM, but this requires a 
more detailed theoretical analysis (EJ, which takes into account the influence of subsurface 
electron scattering on the tunneling current. 

The physical picture of an electron tunneling through a classically forbidden region is that 
the electron flow emerging from the barrier is defined by the matching of the wave functions 
of carriers incident on the barrier and those that are transmitted. For a three dimensional 
STM geometry the wave functions for electrons transmitted through the vacuum region 
radically differs from the electron wave functions in an isolated sample and they describe 
the electron propagation into the bulk from a small region on the surface below the STM tip. 
In contrast, the theory of Ref. [H and its modifications (see [l|, 0] and references therein) 
uses unperturbed wave functions of the surface Bloch states. Changes in the wave functions 
of transmitted electrons due to scattering by subsurface objects provide the information 
about such scattering in the STM conductance. 

In Ref. j6[ it was proposed to introduce in the theory of STM the model by Kulik et al. 
[7j]. In this model a three dimensional STM tip is replaced by an inhomogeneous barrier in 
an otherwise nonconducting interface that separates the two conductors. In Ref. [3] it was 
shown that under assumption of small transparency of the tunnel barrier the wave function 
(and thus the current- volt age characteristics) can be found analytically for an arbitrary size 
of the tunnel area. The results in [fij for the conductance of the tunnel point-contact were 
generalized to an arbitrary Fermi surface for the charge carriers in Ref. In a series of 
papers the model 0] has been expanded to describe oscillations of the STM conductance 
resulting from electron scattering by subsurface defects (for reviews see [9l. [To|| ) . 

Scanning tunneling microscopes have been widely used for the study of various small- 



sized objects: islands, thin films deposited on bulk substrates, etc (ll|-(l8j. First, a discrete 



periodic spatial variation of the STM current originating from the quantization of electron 
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states was observed in the quantum wedge: a nanoscale flat-top Pb island on a stepped 



Si(lll) surface Later these authors showed that that the lattice structure of an interface 



buried under a film of Pb, whose thickness can be as many as 10 times the Fermi wavelength, 



can be clearly imaged with STM [12| . They concluded that the key to the transparency of a 
metal lies in a highly anisotropic motion of the electrons and the strong quantization of their 
transverse wave function components. In the paper 13j the electronic states of thin Ag films 
grown on GaAs(llO) surfaces was investigated by STM with single-layer thickness resolution, 
and the quantum- well states arising from the confinement geometry of the Ag films have been 
identified. Quantum size effects, manifested in the formation of new electronic bound states, 
were investigated by STM on thin Pb islands of varying heights on the Si(lll)-(7x7) surface 



in Ref.|14||. In experiments [15| it was demonstrated that scanning tunneling microscopy and 
spectroscopy of epitaxial Pb islands on Si(lll) reveal adiabatic lateral modulation of the 
energy spectra of the quantum well, providing remote electronic images of the subsurface 
reflection phase. In Ref. [l6j a step structure at the buried Pb on Si(lll) 6x6-Au interface 
was determined by utilizing the presence of quantum well states. It was demonstrated 
that the spatial step positions as well as the step heights can be extracted nondestructively 
and with atomic layer precision by STM. Vertical Friedel oscillations in interface-induced 
surface charge modulations of Pb islands of a few atomic layers on the incommensurate 



Si(lll)-Pb surface have been observed [17|. Thus, detailed experimental results have been 
obtained, but a microscopic theory for STM tunneling spectra on samples of finite size has 
not been reported, which provides the motivation for the present work. Current-voltage 
characteristics for size quantization in planar thin film geometries of metal-insulator-metal 



tunneling junctions have been investigated theoretically in Refs. [19|, |20j. Standing electron 



wave states in thin Pb films have been observed by electron tunneling in early experiments 



by Lutskii et al. [21 



In this paper we present the differential conductance G (V) for small contacts, having a 
radius a smaller than the Fermi wave length Ap = h/pF, where pp is the Fermi momentum. 
The contacts are formed on the surface of a thin metal film and we analyze the voltage 
dependence of I (V) and G(V). We focus on the size quantization effects of the electron 
energy spectrum in the film on G (V). 

The organization of this paper is as follows. The model that we use to describe the contact, 
and the method for obtaining a solution of the three-dimensional Schrodinger equation 
asymptotic in the small radius of the contact, are described in Sec. II. In Sec. Ill the 
current-voltage characteristics and the differential conductance are found on the basis of a 
calculation of the probability current density through the contact. Sec. IV presents a physical 
interpretation of the results obtained. In Sec. V we conclude by discussing the possibilities 
for exploiting these theoretical results for interpretation of electron energy spectroscopy in 
thin films by STM. In the appendixes we solve the Schrodinger equation for the tunnel point 
contact in framework of our model (App.I) and for a point contact without barrier (App.II) 
and find the wave functions for electrons transmitted through the contact. These solutions 
are used in Sec. Ill for the calculation of current. 



II. MODEL AND ELECTRON WAVE FUNCTION OF THE SYSTEM 

The model that we consider is illustrated in Fig. 1. Electrons can tunnel through an 
orifice centered at the point r = in an infinitely thin insulating interface at z = from a 
conducting half-space (the tip) into a conducting sheet of thickness d (Fig. 1(b)). The radius 
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FIG. 1: Schematic representation of a STM experiment on a thin metal film (a) and the model that 
we employ to represent the contact between a bulk conductor (tip) and a metallic film (b). The 
dashed picture of the tip in (a) illustrates a metallic point-contact (STM tip touches the surface). 
Electron trajectories in (b) are shown schematically. 



a of the contact and the thickness d of the film are assumed to be much smaller than the 
shortest mean free path, i.e we consider a purely ballistic problem. The wave function ip 
satisfies the Schrodinger equation 

9m* 

V 2 V (r) + -p- [s - U (r)] i> (r) = 0. (1) 

In Eq.([T]) m* and e are electron effective mass and energy, respectively. The inhomogeneous 
potential barrier in the plane z = we describe by the function U (r) = U f (p) 5 (z), where 
p = (x, y) is a two dimensional position vector in the plane and / _1 (p) = 9 (a — p), with 
(x) the Heaviside step function. For such model the wave function ip (r) satisfies the 
following boundary conditions at the interface z = and at the metal sheet surface z = d 

V>(p,+0)=V(p,-0), (2) 
2m* 1 1 

€ (P, +o) - V4 (p, -o) = (P) V (P, o) , (3) 

^(p,d) = 0. (4) 

Eqs. (JTJ)— (SD can be solved in the limit of a small contact, ka 1 (k = \ / 2m*e/H is the 
absolute value of the electron wave vector k). In the zeroth approximation in the contact 
diameter the solutions of Eq.([T]) for z ^ are independent and satisfy the zero boundary 
condition if>(p, 0) = at the impenetrable interface at z — 0. The quantum states in the 
conducting half-space (z < 0) (the tip) are defined by the three components of the electron 
wave vector k = (ky, k z ), with ky a two-dimensional vector parallel to the interface. In the 
metal film (0 < z < d) the quantum states are characterized by a two-dimensional vector 
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FIG. 2: Space distribution of the square modulus of the wave function for electrons injected by 
a STM tip into a metal sheet of thickness d = 15A, where A = 2n/k is the electron wave length, 
A = A/2vr. 

k perpendicular to the z— axis and by the discrete quantum number n (n = 1,2, ...) which 
results from the finite size of the conductor in the z direction. The energy eigenvalues and 
eigenfunctions for the two disconnected conductors are given by 

h * { k l + k i) = 

2m* 2m* ' 

V> 0) = 2ie i ^ p sin k z z, z < 0; (6) 



and 



^ (r) = -2ie lKp sinfc 2n z, < z < d, (8) 

where fc 2n = 7m/ d. In Eqs. ([6]) and (jSJ) we use a wave function normalization with unit 
amplitude of the wave incident to the interface. 

The partial wave for the first order approximation tp 1 (r) in the small parameter k<l, 
which describes the transition of electrons from one to the other conductor, is given in 
the appendices. Appendix I, Eqs. (AI.5) and (AI.6), gives the solution for a tunnel point 
contact, having a potential barrier of small transparency t = kh 2 /m*Uo <C 1 at the orifice 
in the plane z = 0. Appendix II, Eqs. (AII.6) - (All. 8), gives the solutions for a contact 
without barrier. Figure 2 illustrates the spacial variation of the square modulus of the wave 
function for electrons transmitted through the contact into the film. 



III. CURRENT- VOLTAGE CHARACTERISTIC AND CONDUCTANCE OF A 
POINT CONTACT. 

As has been shown in Ref. 22| for a ballistic point contact of small radius a, with a much 
smaller than the electron mean free path I, the electrical potential V(r) drops over a distance 
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FIG. 3: Illustration of the occupied energy states at zero temperature in the two conductors for 
both signs of the applied bias eV: (a) eV > 0, (b) eV < 0. 



r ~ a from the contact, and in the limit a — > the potential V(r) can be approximated 
by a step function V Q(— z). In this approximation, for the calculation of the electrical 
current we can take the electron distribution functions f^> at z ^ as the Fermi functions 
f F with energies shifted by the applied bias eV (e is the negative electron charge), = 
fp (e — eV 0(— z)) . Figure 3 illustrates the occupied energy states in the two conductors for 
both signs of the applied bias eV. At eV > the electrons flow from the bulk conductor 
(the tip) into the film and, vice-versa, at eV < they flow from the film into the massive 
conductor. The total current through the area of the contact can be found by integration 
over the flux in both directions 



I(V) 



-i poo °° 

/OD 
dkJ< + )/ F (e-eV0(l-/ F (e)), 
-oo 



2nd 
2 



(9) 



(27T) 

In Eq. we integrate over the wave vector k in the semi-infinite conductor for the current 
in the negative direction (first term), and integrate over the two-dimensional wave vector 
K and sum over the discrete quantum number n for the opposite direction of the current 
(second term). 

For simplicity we will take the temperature to be zero. In this case the electric current is 
defined by electrons passing the contact in one direction only, depending on the sign of the 
applied bias. The flux integrated over the area of the contact is calculated in the usual 

way 



h 



2n 



m 



dpp 



>Im 



Pi 



9 , 



(10) 



=±0 
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where p = (pcos0, psin0) . The wave function ij) x (p,z) should be taken as the wave trans- 
mitted through the contact, given by equations (AI.5) and (AII.6) with k = k z for electron 
flux from the tip to the sheet, J^ + \ and by equations (AI.6) and (AIL 7) with Ik = k zn 
[n = 1,2...) for fluxes J^~> in the opposite direction. The energy shift eV in the region z < 
should be taken into account, which for our choice of the reference point of energy (see 
Fig. E]) implies that the absolute value of the electron wave vector in the half-space z < is 
given by k = ^/2m* (e — eV)/h. 

For the tunnel point contact (tpc) the flux can be expressed in terms of the wave function 
in the contact plane (AI.l), and we obtain, 



/(+) „ 

°tpc — 



7r 4 |e| a 4 H 5 k 2 cos 2 $ 
l2m* 3 d 3 U$ 



N{N +1) {2N + 1) 



(11) 



and 



f\ 7i \e\ a 

°tpc — 



4fc5pp 

It l\j tt „„- 



(12) 



Here, $ is the angle between the vector k and the z axis, and N (k) = [kd/ir] with [x] the 
integer part of x. 

For a metallic point contact (mpc) without barrier the expressions for the flux Jmpc are 
written by means of Eqs. (AII.4), (AII.9)-(AII.ll), 



° mpc 



7r 2 \e\ Ha 6 k 2 cos 2 ?? 
9m* d 3 



N{N + 1) (2N+1) 



(13) 



and 



/(-) ~ 

° mpc 



8vr |e| ha 6 k 3 k 2 zn 



9m* 



Substituting Eqs. (ITTT)- (|T^1) into the general expression 
acteristic of the system 



(14) 

we find the current-voltage char- 



I(V) 



(k F dy 



m^_2m^r )s ^ [k) 



F 



eV>0 



(15) 



and 



I{V) 



2 

3 L 




S 5 [k F 



S 2 k F 



(k-pd) 3 
S3 (k F J ~ S 3 (k F 

IT 

kpd 

3 £p 3 
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eV\ 
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(16) 
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FIG. 4: Dependence of the total current, I (V), on the applied bias over the point contact for two 
choices of the thicknesses of the metal film. The constant Iq is given by Eq. (|17p or Eq. f)18|) . 



where e F — fr k F /2m* is the Fermi energy, 



Io. 



tpc 



|e] vr 2 a 4 fr 5 fc F 
12m* 3 [/ 2 ' 



'0,mpc 



eha 6 k F 



9m* ' 

S m (fc) is a finite sum of m — th powers of integer numbers, 



(17) 
(18) 



N(k) 



n 



(19) 



n=l 



Note that S m (k) = H_ m (N) , where H m (n) are generalized harmonic numbers. The cur- 
rent is plotted in Fig. H] as a function of bias voltage for two choices of the film thickness. 
Differentiating Eqs. ( I15p and f|T6|) with respect to voltage we obtain the differential conduc- 
tance G (V) = dl/dV for a point contact with radius a <C X F , 



G (V) =G 1 



k F „ f~ \ 1 



fcp 



fcp 



el/ > 0; 



(20) 



Is 
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|eV|' 






e F _ 



4 / 7T 

3 



S 5 ( k F J - S 5 (fc F ) 



S 2 (k F )+A 
k F 



3h 



l + 4i 



|eV] 7T 
eVl 



£p 



S 3 (k F ) - S 3 (A; F ) 
2" 



+ 



S 2 ( A; F 



(21) 



< 0. 
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In the limit eV — > the zero-bias conductance taken from both sides coincides, as it should, 



G (0) = G X S 2 (hp) = ^N F (N F + 1) (2N F + 1) 

6 



(22) 



where Np = N (kp), and G\ is the conductance of the contact between the bulk conductor 
(the tip) and a thin film that has only a single energy level available below Ep for the motion 
along z, 

S7T 3 

G 1 = G (0)— —3. (23) 



(kpdy 

Go (0) is the conductance of a contact between two conducting unbound half-spaces. For a 
tunnel point contact this is given by 0, 0] 



Go,tpc (0) 



kph 2 \ 2 e 2 {kpa) L 



m*Uo 



36nh 



(24) 



and for a metallic point contact we have [23 



a 



0,mpc 



(0) 



8e 2 (fc F a) ( 
27n 3 h 



(25) 



For d — ?• 00 Eqs. (120]) and (|2ip transform into the known voltage de pen dence of the conduc- 



tance for a point contact between unbound conducting half-spaces [25 



Go {V) = Go (0) 



1 + 



\eV\ 1 f\eV\ 



Ep 



(26) 



The dependence of the differential conductance G (V) for both signs of applied voltage 
is illustrated in Fig. [5J For comparison the dependence G (V) /G (0) from Eq. (j2"B"]) is also 
shown. 



IV. DISCUSSION 

Thus, in the framework of the model illustrated in Fig. [T]we have obtained the current- 
voltage characteristic and the differential conductance for a contact on the surface of a thin 
metal film. Under the assumption that the contact radius a is much smaller than the Fermi 
wave length Xp we found asymptotically exact formulas for the dependence of the total 
current I (V) (Eqs. (JTSj) and (jlBjl ) and the contact conductance G (V) (Eqs. ([20]), on 
the applied voltage. In the limit of zero temperature and neglecting scattering processes 
we have demonstrated that the I ( V) dependence has kinks and G (V) undergoes jumps at 
the same values of applied bias eV (see FigJH and Figj5]) These events result from the size 
quantization of the electron spectrum in the film. 

The results obtained show that even in Ohm's-law approximation ( l22i) . eV — » 0, the 
conductance G (V) is not simply proportional to the electron density of states (DOS) in the 
isolated film, 

( \ m * NF (07\ 

p f {£) = w (27) 
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FIG. 5: Dependence of the normalized differential conductance, G (V) /Go (0), on the applied bias 
over the point contact for two choices of the thicknesses of the metal film. The voltage dependence 
for a point contact between two semi-infinite bulk conductors is shown for comparison (short-dashed 
curve) . 

It is remarkable that the dependence of the conductance G (0) ( 1221) on the number of quan- 
tum levels N-p is the same for, both, tunnel and metallic point contacts. This fact shows that 
such dependence is not sensitive to the model taken for the potential barrier, and that it is 
the result of the point-contact geometry. Recently, the relationship between the differential 
conductance and the local density of states has been studied in a tight-binding approxima- 
tion for tunnel junctions, where the junction geometry can be varied between the limiting 



cases of a point-contact and a planar junction 26]. In the framework of a real-space Keldysh 



formalism the authors of Ref. [26] have shown that the differential conductance is not, in 
general, proportional to the sample DOS for planar junctions, although features of the DOS 
may be present. 

From Eqs. ( 12"U|) and ( l2~Tj) it follows that the the conductance is non-symmetric in the 
applied bias. This asymmetry can be explained as follows: Let eV > and electrons 
tunnel from the bulk conductor into the film (Fig. [3^) in which Np subbands of the size 
quantization are partially filled. If the bias eV is smaller than the distance Ae between the 
Fermi level Ep and the bottom of the next (empty) subband En+i = 7f 2 h 2 (Np + l) 2 /2m*d 2 , 
Ae = en+i — £f, the electron can tunnel into any of the Np subbands. At eV = Ae tunneling 
into the (N-p + l)-th subband becomes possible and the conductance G (V) undergoes a 
positive jump. Such jumps are repeated for increasing voltage for all higher subbands. For 
eV < 0, when electrons tunnel from the thin film into bulk metal (Fig. [3b) the situation is 
somewhat different. If the bias \eV\ becomes larger than distance Ae between the bottom 
of the last partially filled subband En = 7r 2 h 2 Np/2m*d 2 and Fermi energy, Ae = Ep — En, 
the contribution of the iVF-th subband to the tunnel current does not depend on the voltage 
because for any \eV\ > Ae all electrons of this subband can tunnel into the bulk states 
of the left conductor. For this reason the differential conductance drops for values of \eV\ 
coinciding with bottoms of subbands of size quantization in the film. The distance between 
neighboring jumps of the conductance on the voltage scale equals the distance between 
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energy levels Aen = £n+i — en = n 2 h 2 (2N F + 1) /2m*cP. For eV < the number of 
conductance jumps is finite and equals the number of discrete levels below Fermi surface 
Np. The asymmetry around V = and the general shape of the jumps in the conductance 
can be recognized in the experiments, see e.g. [l3j]. 111 t ne special case of a 2D electron 
system, which has only one level in the potential well, there is a single negative jump of 
G (V) . Such a jump has been observed in Ref. [27j by STM investigations of the 2D electron 
gas at noble-metal surfaces. For eV > the number of conductance jumps formally is not 
restricted. However, for eV > Ep our approach is no longer applicable and the influence of 
field emission on the tunnel current must be taken into account [28, 29 



The observation of manifestations of the size quantization in the STM conductance re- 
quires a few conditions which must be fulfilled: The distance between the energy levels must 
be large enough and should satisfy the condition Ae n ^> h/r,T, where r is the mean scat- 
tering time of the electrons in the film and T is the temperature. The surfaces of the metal 
film in the region of the contact must be atomically smooth [30]. When a finite lifetime of 
the quantized states becomes relevant, the temperature broadening of the Fermi function, 
or surface imperfections need to be taken into account this will result in a rounding of the 
jumps in the curve G (V) presented in Fig. [5] (Eqs|20j [21]), which was plotted under as- 
sumptions of perfectly specular surfaces, T = 0, and r — )■ oo. With these restrictions taken 
into account the current- voltage curves in Fig. H] give a fair qualitative description of the 
experimental results of Ref. [12j . 

It can be easily seen that the results obtained have a more wide domain of applicability 
than that of a rectangular well for the conducting film. For any model of the potential 
which restricts the electron motion in one direction the differential conductance has a step- 
like dependence on the applied bias with distances between the steps equal to the distances 
between the quantum levels. 



V. CONCLUSION. 

In summary, we have investigated the conductance of ultra small contacts, for which the 
radius is smaller than the Fermi wave length, on top of the surface of a thin film. The dis- 
creteness of the component of the electron momentum transverse to the film surface is taken 
into account, where the distance between the electron energy levels due to the size quantiza- 
tion is assumed to be larger than the temperature. Both, a contact with a potential barrier 
of low transparency, and a contact without barrier have been considered. In framework of 
our model, using a 5-function potential barrier, the current-voltage characteristic / (V) of 
the system and differential conductance G (V) have been obtained. We predict a sawtooth 
dependence of G (V) on the applied bias and show that the distance between neighboring 
jumps is equal to the distance between neighboring energy levels of size quantization, i.e. 
this dependence can be used for spectroscopy of size quantization levels. At eV > the 
jumps in the conductance are positive and correspond to distances between levels above the 
Fermi surface, while G (V) undergoes negative jumps for eV < 0, the distances between 
which are equal to the distances between the levels below the Fermi surface. The predicted 
quantization of the conductance can be observed in STS experiments, and the shape of the 
theoretical curves agrees well with experiments. 
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VI. APPENDIX I: ELECTRON TUNNELING BETWEEN THE TIP AND THE 
THIN FILM. 



We search a solution to Eq. (pQ) at V = in the form of a sum ip = ip + ip 1 for the 
incident and backscattered waves, and ip — for the transmitted wave. Here, ip , as given 
by Eqs. (jBJ), flE}, is the unperturbed wave function that does not depend on the barrier 
amplitude U , while ip 1 ~ 1/U gives the first order correction. Substituting the wave 
function into the boundary conditions ([2]) and fl3]) one should match terms of the same order 
in 1/Uq. As a result the boundary condition (J3J) becomes 

Vi(p,0) = ^6(a-p), (AI.l) 

where k = when the wave is incident to the contact from the tip side, and k = k zn when 
the wave arrives at the contact from the sheet. For fca C 1 we have in the plane of the 
contact np <C 1 and we can neglect the exponent in the boundary condition (AI.l). 
The function ip l (p, z) can be represented as a Fourier integral 

oo 

V>i (p, z) = J dn'e~^' p ^ (*', z) . (AI.2) 

— oo 

The Fourier components in (AI.2) should satisfy the zero boundary condition at z = d, but 
are otherwise freely propagating along z, 

V( K ',z) = y(K',0) Smk ' z ( f ~ d \ O^z^d, (AI.3) 

sm k' z d 

z) = 0) exp {-ik' z z) , z^O, (AI.4) 

with k' z = yjk 2 - k' 2 , k = y/2m*e/h. From Eqs. (AI.l), (AI.2) it follows that 



OO 



. , , 1 /',,•-'„,/ v ikh a J\ (k'cl) 
*< K ' 0) = Wf J = - 2«m>U V - (AU) 



-OO 



Substituting this into Eq. (AI.2) we find the wave functions for the electrons transmitted 
through the contact as, 



oo 



, . , ikh a /"„,,,, , , , sin ki ( d — z ) 

^ (p, z) = — - / dn'J (k' P ) J, ( K 'a) . z \ \ < z < d, (AI.5) 

o 

oo 

ik^ 2 a f 

ip-, (p, z) = — — / dn'Jo (k'o) Ji in' a) exp (—ik'z) , z < (AI.6) 

m*[/ J 



o 



where J„ fx) is the Bessel function of the first kind. 
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VII. APPENDIX II: METALLIC POINT CONTACT BETWEEN STM TIP AND 
METAL FILM. 



Here we consider a point contact without potential barrier in the plane of the interface. 
When the contact radius is small, ka <C 1, we can use perturbation theory for the electron 
wave function in the limit a — > 0. In zeroth approximation the wave functions are given by 
Eqs. (jBJ) and (jSJ). The first order correction, ip 1 (p, 0), to the wave function in the plane of 
the contact can be found by the method proposed in [23j |. For distances r < A from the 
contact we can neglect the second term in the Schrodinger equation (P} and it reduces to 
the Laplace equation. We express the wave function in coordinates of an oblate ellipsoid 
of revolution (a,T,(p), with a ^ and — 1 ^ r ^ 1. As a consequence of the cylindrical 
symmetry of the problem the wave function if> 1 (a, r) does not depend on <p. The interface 
corresponds to r = and the plane of the orifice is at o = 0. In these coordinates we obtain 
the equation 



d_ 



(1 



a 



da 



d 



(1 - r 2 ) 



dip! 
dr 



0. 



with the boundary condition at the interface 

if> 1 (a>0,T = 0) = 0. 
The solution of the boundary problem (AII.l), (AII.2) is 

(a, t) — t [c\G + C2 (1 + a arctan a)] 



(AII.l) 



(AII.2) 



(AII.3) 



where C\ and are constants. For a = Eq. (AII.3) gives the function ip 1 (p, z) in the plane 
of the contact z = 0, p ^ a 



ipt (p,0) = c 2 Wl-^. 



(AII.4) 



As in Appendix I, we express ip l (p, z) as a Fourier integral and, using the Eq. (AII.4), we 
find for the Fourier components, 



(k',0) 



1 



(27T) 5 



dpe iK ' p ip 1 (p, 0) = c 2 a 



3i (k a) 



K 



(AII.5) 



where j\ (x) is the spherical Bessel function of the first kind. Substituting Eq. (AII.5) into 
Eq. (AI.2) and using Eqs. (AI.3), (AI.4) we obtain 



c 2 a 
27 



dn'J (n'p) ji (k'ci] 



sin k' z (d 



sin k' z d 



< z €d, 



(AII.6) 



and 



c 2 a 
2^ 



dn'Jo (k'p) ji (n'a) e 



The constant c 2 must be found from the boundary condition 
case takes the form 

(9?Mp,+0) «9^(p,-0) 



dz 



dz 



2ik = 0. 



z < 0. 



(AIL 7) 



at Uq = 0, which for this 
(AILS) 
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The meaning of the symbol k is explained below Eq.(AI.l). Differentiating Eqs. (All. 6) and 
(AII.7) with respect to z and calculating the integrals in the limit of small a we find, 



dz 



z=+0 



c 2 a 
~2n 



IT 



n 3 a 



±T l-Trz+iTTrzW + WN+l) 



2a 2 18d 3 



dz 



c 2 a 
2t7 



IT 
2^ 



.irk 3 a 
i 

9 



(AII.9) 



(AII.IO) 



where A" = [kd/ir] with [x] the integer part of x. Substituting Eqs. (AII.9) and (AII.IO) into 
(AIL 8) in leading approximation in a, in which only first terms in brackets (proportional to 
1/a 2 ) should be taken into account, we find for the unknown constant, 



c 2 ~ 2ika. 



(AII.ll) 
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